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Abstract 

We calculate the spectral function A p of the p meson in nuclear 
matter. The calculation is performed in the low density approxima- 
tion, where the in-medium self energy S mef ^ is completely determined 
by the vacuum pN forward scattering amplitude. This amplitude is 
derived from a relativistic resonance model. In comparison to previous 
non-relativistic calculations we find a much weaker momentum depen- 
dence of £ me d, especially in the transverse channel. Special attention 
is paid to uncertainties in the model. Thus, we compare the impact 
of different coupling schemes for the RNp interaction on the results 
and discuss various resonance parameter sets. 
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1 Introduction 



The question of how vector meson properties change at finite baryon density 
or temperature has triggered a lot of theoretical and experimental work in 
the last decade. Of particular interest are the in-medium properties of p 
and ax meson, which are chiral partners. One of the most evident signa- 
tures of the spontaneous breaking of chiral symmetry in vacuum is the large 
mass difference between these mesons. The claim is that the in-medium 
mass spectra of p and a% meson can be related to the partial restauration 
of chiral symmetry. If a chirally symmetric phase is reached at high den- 
sity/temperature, a degeneracy of their spectral functions is expected [fl], £3]. 
Whereas only a few studies exist on the in-medium properties of the a\ me- 
son H, U, numerous theoretical works have been performed for the p meson 
0, |, % 1, g 0, 0, 0, pi, §0. 

It has been argued by several groups that the mass of the p meson is 
related to the expectation value of the scalar quark condensate. The latter is 
supposed to decrease at finite nuclear density due to the restauration of chiral 
symmetry, thereby inducing a mass shift on the p meson. The most famous 
of these models is that of Brown and Rho ||. It predicts a downward shift 
for the p meson mass of about 15% — 20% at normal nuclear matter density. 
Similar results have been found in |6j using the QCD sum rule approach 
under the assumption that the width of the p meson remains small in the 
nuclear medium (cf. for a thorough discussion of that point). 

In other approaches the properties of the p meson in nuclear matter are 
calculated by taking into account conventional nuclear many body effects. At 
low nuclear densities p^ the self energy T, med can be calculated in the p^T 
approximation (low density approximation) and is completely determined by 
the knowledge of the vacuum p N forward scattering amplitude T. Since this 
quantity is not directly accessible by experiment, one needs to construct a 
model which incorporates experimental constraints as far as possible. This 
was first pointed out in M. Information on pN scattering comes mostly 
from tcN — > ttttN processes. One of the most successful descriptions of these 



reactions has been achieved by Manley et al within an isobar model [18], 

A resonance model for the pN scattering amplitude allows the incor- 
poration of experimental constraints in a straightforward fashion. Such a 
strategy was first suggested in |9[], where only p— wave resonances were taken 
into account. However, as shown in our previous paper flOU , the _Di 3 (1520) 
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resonance plays a dominant role in pN scattering and can therefore not be 
neglected in a complete analysis. In both models the resonance parameters 
were taken from the PDG J2TJ and the calculation was performed within a 
non-relativistic framework. An alternative approach has been chosen in [|TT| 



where the baryonic resonances are not introduced as excited three-quark 
states, but are generated dynamically in a coupled channel approach from 
meson-nucleon and meson-A rescattering. In [|l^] the pN scattering ampli- 
tude is derived from an effective chiral Lagrangian theory, neglecting, how- 
ever, resonant contributions. 

Going beyond the low density approximation, various groups have con- 
sidered the change in the 2 7r decay of the p as induced by the modified tt 

H31 



0, 0, |I7J. In |T§ [T7 



spectrum in the nuclear medium, see e.g. 
resonance contributions have also been included. 

Qualitatively, all these models predict a strong broadening of the p me- 
son. The mass, however, remains nearly unchanged. This finding is in good 
agreement with QCD sum rules, see e.g. [0, [L^]. 

So far experimental signals indicating a change of the p spectrum in the 

, |23| and the 



nuclear medium, have been obtained by the CERES |2T], 

In these experiments the dilepton spectra display 



HELIOS collaboration M 



a downward- shift of strength in the region of the p meson as compared to 
theoretical predictions. Agreement between theory and experiment can be 
reached by introducing medium modifications for the p meson [25]. On the 
other hand, in [2£| the claim is made, that the CERES data do not necessarily 



hint to a medium modification of the p meson, but can also be explained with 
its vacuum properties. 

Certainly, in order to gain a more quantitative understanding of the p 
meson properties in nuclear matter, further theoretical studies, especially on 
the momentum dependence of the spectral function, and improved experi- 
mental statistics, are mandatory. The sensitivity on details of the spectral 
function is even enhanced in photo-nuclear reactions, as was demonstrated 
in 0. 

[TO 



m Ira 



have treated the res- 
As we will show, this in- 



Up to now all existing calculations 
onance contribution in a non-relativistic manner, 
troduces substantial uncertainties concerning the momentum dependence of 
the p spectral function. We therefore do not only lift the non-relativistic 
approximation, but also study the dependence of the results on various cou- 
pling schemes. Furthermore, we analyse the influence of different resonance 
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parameter sets on the results. We thus examine up to which extent safe con- 
clusions about the propagation of p mesons in nuclear matter can be drawn 
from a resonance model. 

The paper is organized as follows. In Sect. [|we will present an overview 
of the formalism employed to describe the spectral function. Furthermore, 
an arbitrariness in the non-relativistic approach is exposed and thus the need 
for a relativistic model is motivated. The basic ingredients of the relativistic 
calculation will be introduced in Sect. |3|. A problem related with the de- 
scription of the resonance propagator and its solution is presented in Sect. 
|3. 1| . In Sect. (| we analyze the quality and properties of different resonance 
parameter sets and their influence on the results. In Sect. [| the results 
of the relativistic calculation are presented and then compared to the non- 
relativistic calculation. The impact of various coupling schemes for the vertex 
RNp is studied in Sect. We summarize our results in Sect. |7]. 



2 Overview of the Model 

In this section we explain the formalism on which our calculation of the in- 
medium spectral function of the p meson is based and introduce some basic 
notation. We then review the non-relativistic scheme which has been em- 
ployed in all works up to now || [lO|, [Tj^ |17| and use some simple kinematical 



considerations to point out why a fully relativistic framework is needed. 



2.1 Definition of the Basic Quantities 

The spectral function is defined as the imaginary part of the propagator. Due 
to Lorentz invariance the vacuum spectral function of the p meson depends 
only on one kinematic variable - its invariant mass - and both transverse 
and longitudinal modes have the same spectral distribution. The presence 
of nuclear matter breaks Lorentz invariance. Consequently transverse and 
longitudinal propagation modes of the p meson receive different in-medium 
modification. Also, the corresponding spectral functions depend both on the 
energy uj of the p meson and its three-momentum q: 

A W(U, q) = 1 ImST/ >^ 

p 1 ' 4J 7r(cj 2 -q 2 -m 2 + Re£ T / L (u;,q)) 2 + ImE T / L (w,q) 2 

(1) 
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m p is the pole mass of the p meson. £ T/,L is the transverse or longitudinal 
part of the self energy of the p meson and can be decomposed into a vacuum 
and an in-medium part: 

S^(c,q) = S_(g 2 ) + S^(c,q) . (2) 

The imaginary part of the former is to lowest order in the coupling con- 
stant given by the 2-pion decay width of the p meson [ 12] , |14| : 



Im E vac (q ) = y/q T n7T (q ) , (3) 

and 

2 / \ 3 

(4) 




Here q is the squared invariant mass of the p meson, T its decay width on 
the pole mass and p mp (p) denotes the 3-momentum of the pions measured 
in the rest frame of a decaying p meson with mass m p (y/q 2 )- For invariant 
masses of the p meson up to around 1 GeV - which we consider in this paper 
- the real part of £„ ac (g 2 ) shows only small variations and can be neglected 

H- 

At finite nucleon density p N the in-medium self energy is given by 
the following expression: 

^ T / L t \ — f ^PN rpT/L, 



In our notation q = (a?, q) and = (E N , p N ) are the 4-vectors of p meson 
and nucleon, respectively. The integration is restricted to the volume of 
the Fermi volume, denoted here by Q. By performing the integration in 
the rest frame of nuclear matter Q is simply the Fermi sphere, and the pn 
integration reaches from to the Fermi momentum pp. T t J t is the pN 
forward scattering amplitude, as discussed below. Note that T, me d(uJ, q) is a 
Lorentz invariant quantity, since both the forward scattering amplitude T^( L 
and the integrations measure ^F^Em are Lorentz scalars. 

At sufficiently low nuclear densities it is advantageous to expand the self 
energy in terms of the Fermi momentum. This is most easily done in the rest 
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frame of nuclear matter. The 4- momentum of the nucleon simply becomes 



Pn = (tun, 0). One then recovers the pn T or low density approximation [pg 



s med(9) = ^—PNTH L (q,p N ) . (6) 

The forward scattering amplitude T^J t L is obtained as a sum over all in- 
cluded resonances. For a single resonance the scattering amplitude is depicted 
in Fig. |I|. The corresponding contribution to the transverse and longitudinal 
channel T T I L is obtained from the tensor T^ u (q,p^) upon contracting with 
the projection operators P^J L [0, [T2|, H| : 



t t = _ T „u pT 
2 ^ 



(7) 



The general form of T^ u (q, p^) reads: 

r f 2 TrV?(q,p N )Di(k)VZ(q,p N )W + m N ) for J = § , 

^rrVT^P^^^Vf^p^^ + m^) for J = f , 

v 2 2 f 2 

(8) 

where k = p^ + q is the four- momentum of the intermediate resonance. The 
coupling strength at the RNp vertex is given by /. A detailed discussion of 
both the relativistic vertex functions Vl(q,pN) and V^{q,pN) as well as the 

2 2 

relativistic propagators of spin-i and spin-| resonances Di(k) and Df(k) 

2 2 

will follow in Sect. [3]. 

The coupling constant / is obtained from a fit to the measured partial 
decay width R — > N p. In Sect. (| we discuss in detail from which analysis 
the resonance parameters are taken. Since many of the resonances which are 
considered in this work are below or relatively close to the nominal threshold 
for this decay of yfs = m N + m p ~ 1.7 GeV, we integrate over the vacuum 
spectral function of the p meson (see Fig. ||]). The decay width then reads: 

\/s~m N 

^n p (s) = 7^—71^- I dm2mF(s) 2 
2j R + 1 8ns J 

x ^ C (m 2 ) (2 |-M^p| 2 + \M L RN /) qcm • (9) 
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The spin of the resonance is denoted by jn } Ir is an isospin factor. A v p ac 
is the vacuum spectral function of the p meson. By q cm we denote the 3- 
momentum of the p meson in the rest frame of the decaying resonance. F(s) 
is a formfactor that will be explained in Eq. [21]. The squared matrix element 
for the resonance decay |-M^^| 2 reads in terms of the vertex functions V 
(here for spin-^ resonances): 

\M T R f p \ 2 = -f RNp 5£ E n <Pn) %(pn) V| ««(*) . 

(10) 

The spinors u^{k) and u s {j>n) denote resonance and nucleon spinors respec- 
tively. The numerical values for the coupling strength fnN P , which we obtain 
this way, are given in Table |l|. We have defined Tn p for an arbitrary invariant 
mass \fs of the resonance, since we use Eq. |9| also as a parameterization of 
the mass dependence of the decay width. 

The pnT approximation Eq. || only describes two-body processes in the 
nuclear medium. At sufficiently large densities also higher order effects need 
to be taken into account - currently, however, a scale that pins down the 
breakdown of the p^ T approximation is not at hand. A comprehensive 
description of many-body effects is very hard to obtain [[H], [L7|]. Therefore, 
we incorporate only two typical in-medium corrections to the low-density 
approximation: 

• integration over the Fermi sphere of the nucleons, 

• renormalization of £ me d due to short range correlations. 

In the case of p-wave resonances, the short range correlations are incorporated 
via the Landau-Migdal parameters: 

E med (u;,q) = Sr " e s d . (11) 

1 — g'^^- 

Only the P33(1232) is renormalized, leading to a suppression of its contribu- 
tion. For the Landau-Migdal parameter we take g' = 0.5. For higher p-wave 
resonances this parameter is not known. The actual calculation shows very 
little effect of g' on the results. For s-wave resonances there is no informa- 
tion available on how to include short range correlations; thus we do not 
renormalize the contribution from these resonances. 
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2.2 Critical Review of the Non-Relativistic Model 



We quickly review here the qualitative features of the results of a non- 
relativistic approach as found in our previous work ]l(J . At low momenta the 
results are completely determined by the properties of the excited resonances. 
Their quantum numbers and experimentally observed coupling strengths to 
the pN channel provide the necessary information. Close to threshold the 
resonance contribution scales like q 2/ , where q is the 3-momentum of the 
p meson. Consequently, at small values of q s-wave resonances (P = — 1) 
yield finite contributions, whereas the contribution from p-wave resonances 
(P = +1) is negligible. Quantitatively, the results of our explicit calculation 
of ^ m ed [TO revea l that the propagation of a p meson in nuclear matter is 
dominated by the s-wave resonances _Di 3 (1520) and -D 33 (1700) as a direct 



consequence of their large coupling to the p N channel as given in p0[ . The 



spectral function receives a rich structure: a new peak arises at an invariant 
mass of about 0.5 GeV from the excitation of the _Di 3 (1520) and the peak at 
the pole mass of the p meson is broadened by the .D 33 (1700). 

Whereas at low momenta the main uncertainty of the results comes from 
the values of the coupling constants, at high momenta an additional model 
dependence is introduced by the momentum dependence of the self energy. In 
a non-relativistic calculation p-wave resonances do not couple to longitudinal 
p mesons. Furthermore, A T and A L display a completely different behaviour 
at high momenta. For A L the effects from resonance excitation are rather 
small and the spectral function is strongly peaked at the p mass. A T , in 
contrast, is completely smeared out in this kinematic regime. 

We will demonstrate now that the momentum dependence of the self 
energy can not be reliably determined within a non-relativistic calculation 
and we argue that therefore a relativistic framework is imperative in order 
to put the results on a more solid basis. 

The forward scattering amplitude is proportional to f 2 V 2 , see Eq. |[ 
From Eqs. |] and |10| it is clear that the vertex function V enters also into 
the determination of the the coupling constant /. It thus appears twice in 
the forward scattering amplitude. An ambiguity now arises in the choice of 
the reference frame in which the non-relativistic matrix element is evaluated. 
In all previous works 0, [H], [16], [T?| the following choice was made: the 
coupling constant was obtained by evaluating the matrix element in the cm- 
frame (rest frame of the resonance) whereas the self energy was calculated 
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in the lab- frame (rest frame of nuclear matter). This is a suggestive choice 
since the decay width is usually defined in the rest frame of the decaying 
particle whereas a nuclear matter calculation is preferably performed in the 
rest frame of nuclear matter. The non-relativistic reduction itself does not 
uniquely determine the reference frame, as it assumes that both the resonance 
and nucleon momentum are small compared to their masses. This condition 
is equally well met in the cm-frame and in the lab-frame. 

As shown in 0, |1(J , for a resonance of positive parity both the resonance 



decay width and the forward scattering amplitude are proportional to q 2 . For 
a given invariant energy y/s the three-momenta in lab-frame and cm-frame 
are easily related: 

y/s 

<ilab = <lcm • (12) 

m N 

For the decay width T one finds r ~ f 2 q 2 . Since the coupling constant / 
is determined from the value of F for an on-shell resonance (see Eq. |J), one 
has (in obvious notation): 

T ~ flab rfab = fern Qcm > ( 13 ) 

thus leading to the relation 

2 

flab - fem ~~T ■ (I 4 ) 

m R 

Similarly, the forward scattering amplitude is proportional to / 2 q 2 . The 
most relevant contribution from each resonance to the self energy comes 
from the kinematical region yfs ~ tur, where 

Ti a b = T cm = T , (15) 

as long as lab-frame or cm-frame kinematics are used consistently. However, 
applying mixed kinematics as in || [H| [16], |17 |, where the coupling constant 



is determined within cm-kinematics and the forward scattering amplitude in 
lab-frame kinematics, produces 

2 

Tmix = fem ^Lfe = T J- ■ (16) 

m N 
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Typically, the resonance masses vary between 1.5 and 2 GeV. Thus for p-wave 
resonances the ambiguity in the non-relativistic reduction leads to uncertain- 
ties which can be as large as a factor of four. For s-wave resonances the effects 
are less pronounced. Here the transverse and the longitudinal channel are 
proportional to the energy uo 2 and the invariant mass q 2 of the p meson, 
respectively. Thus only the transverse channel will be affected by the ambi- 
guity of the reference frame. Noticeable uncertainties for the s-wave part of 
Timed are expected for three-momenta comparable to the invariant mass of 
the p. Altogether, this is clearly not a satisfying situation which can only be 
resolved within a relativistic approach. 



3 Relativistic Framework 



In this Section we set up the relativistic framework for the calculation of S me ^. 
We will present the Lagrangians used to describe the RNp interaction. A 



discussion of the resonance propagator will be deferred to Sect. |Q 



The vertex functions and V3" are derived from the following relativistic 

2 



and gauge invariant Lagrangians: 



C 



hit 





for 


J w = 


3- 
2 




for 


r = 


3+ 
2 




for 


j n = 


1- 
2 


l^UL^a^iliF^ + h.c. 


for 


r = 


1+ 
2 



(17) 



Here ip^ 



denotes the resonance spinor and ip the nucleon spinor, a^ v 



\ Yl^yl v \ an d F^ v = p u — d u p M . We note in passing that from these 
Lagrangians the non-relativistic expressions from our previous publication 
|lOf| are readily obtained. We require gauge invariant interactions in order to 
ensure physical results in the vicinity of the photon point q 2 = 0. 

The coupling of spin-1 particles to the baryon resonances is not uniquely 
determined by its quantum numbers. In analogy to the commonly used 
couplings for the A iV 7 system, we employ - in addition to the couplings in 
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Eq. [T7| - the following couplings in the p-wave channel (see for comparison 



e.g. 



C lnl = { ' ! (18) 
Isgf.^'f^&'Fr + h.c. 

In general, the dimensionless coupling constant J'rnp is different in all 
three cases. For resonances with negative parity, the corresponding La- 
grangians Ci nt can be obtained by dropping the 7 5 -matrix. In Sect. ||we will 
explore the sensitivity of the resonance model on the chosen vertex function. 
For the comparison between non-relativistic and relativistic calculation in 
Sect. |5| the coupling from Eq. [IT] will be employed. 



In view of the difficulties in constructing a relativistic theory for spin-| 
fields, we treat them non-relativistically. In the relevant kinematical region 
only resonances with positive parity couple sufficiently to the pN system. 
We can therefore restrict ourselves to the discussion of p-wave channels. In 
agreement with previous studies ||, [l(J the following Lagrangian is used: 



C int = ^xiR ij F ijX + h.c. . (19) 
m p 

Here Xr an d X are resonance and nucleon spinor. R l i denotes the spin- 
| transition operator ||. We will come back to the treatment of spin-| 
particles in Sect. |5] and Sect. ||. 

For the isospin coupling we take the standard form: 



Xr Ok PkX for I =\ 
Xr S k PkX for 1=1 



(20) 



Xr. and x denote resonance and nucleon isospinor, respectively, S& is the 
isospin- 1 transition operator and the standard Pauli matrix. 

In the actual calculation each vertex is multiplied by a phenomenological 
formfactor [p9~| : 

f M = A , + (s A l m ^ • < 21 > 
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This formfactor has a very different functional dependence on s compared to 
the one previously used in the non-relativistic treatment: 



A 2 



ah) 



\ A 2 + qfU 



(22) 



We reject this formfactor because it is not Lorentz invariant. Furthermore, 
since the resonance excitation is a typical s-channel process, the formfactor 
should have a maximum on the resonance mass-shell and should only be a 
function of the kinematical variable s. Both criteria are obviously not fulfilled 
in Eq. ^ since the three-momentum is a function of the invariant mass q 2 
of the p meson: 



2 

Qlab 



m 



N 



q 



2\2 



4:171% 



q 



(23) 



For the same reason a formfactor depending on the cm-momentum of the p 
meson like 



F{q cm ) 



A 2 



QOcm 



\ A 2 + q! 



2 

cm 



(24) 



is also not satisfactory q C m and qo cm depend both on q . Here qocm 
denotes the p momentum in the cm-frame if the decaying resonance is on its 
mass-shell. 

For the cutoff parameter we take a value of A = 1 GeV for all resonances, 
in accordance with | 29fl . The influence of the new formfactor Eq. |2l| on the 
results will be discussed in Sect. [|. 



3.1 The Resonance Propagator 

In this Section we discuss the propagator of the formed baryon resonance. 
A proper description of this quantity is not a trivial task since it displays 
in general - already in the spin-| sector - a very rich spinor structure (see 
e.q. [pl| ). We show that a commonly used recipe for both the propagators of 
spin- 1 and spin-| particles |32], [33[ leads to unphysical results in the case of 
p N scattering. Then we derive an improved propagator which - albeit being 
only a crude approximation of the exact solution - is nonetheless sufficient 
for our purposes. 
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Unitarity relates the matrix element for the decay of a resonance R into a 
nucleon and a p meson of given polarization [-M^^l 2 to the imaginary part 
of the pN forward scattering amplitude T T I L : 

~ r p r toi (25) 

~ \M T J N L P \ 2 • 

It may be objected that for tree-level calculations unitarity is by no means 
guaranteed [pER . However by considering only the contribution of one res- 
onance in the s-channel, unitarity can be fulfilled provided that the decay 
width and and the resonance propagator are parameterized in a consistent 
way. 

Let us start with the calculation of the squared matrix element for the 
resonance decay |A^^ p | 2 . The standard evaluation of this quantity, that 
contains a trace over the internal degrees of freedom, involves the complete- 
ness relation of the resonance spinor. Taking it naively as 

5>?(fc) =jt + m R (26) 



even for resonances off the mass-shell, i. e. with v k 2 tur, leads to negative 
values for the squared matrix element for positive parity resonances, which 
is obviously wrong. The calculation leading to this results is given in the 
Appendix. Thus, Eq. BB is not appropriate in the off-shell case and needs to 



be modified. The obvious choice is to set [36, 37 



X>?(Jfe) u?(Jfc) = ft+Vk^ . (27) 

r 

Unitarity now enforces that the modified completeness relation Eq. |27| 
also changes the functional dependence of the (imaginary part of the) prop- 
agator. For a spin-| particle one thus obtains 



Di(k) = * , y= (28) 

2 k 2 — m R + iyk 2 T 

instead of 

D? ial (k) = ft + ^fl (29) 

2 k 2 -m 2 R + iVk?T 
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The latter form is often employed in the literature Since - from 

unitarity - it is related to the completeness relation with rriR, it leads to 
negative values for the total pN cross-section 

In the Appendix we give some theoretical backing for Eq. |28|. In particu- 
lar, we show that Eq. |29| violates a fundamental constraint for the propagator 
and should therefore be rejected. We also discuss that Eq. ^ is only correct 
for the imaginary part of Di(k), but not for the real part. In an explicit 
calculation we have checked nowever, that the differences between the exact 
solution and Re Di (k) as taken from Eq. [28| are negligible. 

The propag ator Df(k) of non-interacting spin-| particles is given by the 

2 

well known Rarita-Schwinger propagator [pgfl . The conventional choice of 
adding a finite decay width in the denominator J33], Q leads - as in the case 



of spin- 1 particles - to negative total cross-sections for processes including 



intermediate J 77 



| resonances. As explained in the Appendix we propose 
therefore that the propagator should read: 



Df(k) 



jk+Vk 2 



k 2 



x 
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m 



ai3 



1 



k 2 T 



R 2k a k? 1 k 7^ — k p 
_ T " 7 p 1 ! — 

3 ' 1 3 k 2 3 v/p 



r 



(30) 



Some remarks about this propagator are in order. First we note that 
it is directly proportional to the spin-| projector P3 , see ||35|| . Second, in 

analogy to the case of spin-| resonances this expression is - strictly speaking 
- not correct for the real part, but again we expect the deviations to be 
negligible for our purposes. Finally, Eq. |30| has already been discussed in 
the literature |35|, In [^] the replacement rriR — > \fk? j is introduced 

to insure gauge invariance. This is not necessary in our case as the vector 
mesons are introduced via the gauge invariant field tensor F^ v . On the other 
hand, in |55) it was shown that Eq. [3(] is incorrect in the limit of a small decay 
width T. This observation is in full agreement with our previous statements 
about the real part of Eqs. |28] and |30|. Also, as pointed out in |35], |33|, the 
propagator Eq. |3^ displays a pole at k 2 = 0, which can lead to divergent 
results for scattering amplitudes. However, this kinematical regime is reached 
only in w-channel processes which we do not include in our model. Therefore 
the results of this work should not suffer from these poles. One should also 
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be aware of the fact that the use of gauge- invariant couplings for the spin-| 
field - which introduce additional orders of the resonance momentum - as 
suggested in pi avoids these poles. 



4 Included Resonances And Their Parame- 
ters 

In this Section we list the resonances included in our model for T* me d and 
discuss their parameters such as mass, decay width and branching ratios. 

We include all baryon resonances with a sizeable coupling to the pN 
channel. They are listed in Table |I[ The nucleon is no longer included in 
the model. We found in a relativistic calculation on the T-matrix level that 
the Born terms lead to a strong overestimation of the data for the reactions 
7 N —>■ p N and tt N — > pN, as well as for reactions with an u meson in 
the final state. It is very likely, that on this level important backscattering 
effects are missing, since the nucleon has a large coupling strength to the 
Ntt, N uj and N p channels. Similar effects are known from Compton scat- 
tering, where in a i^-matrix analysis, which includes the rescattering, the 
Born terms yield satisfactory results. In a T-matrix approach on the other 
hand, additional formfactors, which suppress the contribution in backward 



direction, have to be introduced in an ad-hoc manner |R3] . Therefore it does 



not seem appropriate to include the nucleon in our model. For resonances, 
on the other hand, backscattering effects are effectively taken into account 
by their decay width. 



In contrast to our previous publication l lQfl , we take the resonance pa- 
rameters from the analysis of Manley et al [jT^, [19j rather than those from 
the PDG group [pOfl . The PDG parameters represent estimates of resonance 
properties, averaging the available experimental results. They are based on 
a variety of theoretical models, which in general differ in many details, as, 
for example, the included resonances or the treatment of the non-resonant 
background. Therefore, the extracted parameters are strongly model depen- 
dent and it is questionable whether an average of the results from different 
studies provides conclusive information on any specific reaction channel. To 
ensure a realistic description of the p N scattering amplitude, it is much safer 
to stick to the results of a single analysis of tt N — > n n N reactions. In the 
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work of Manley et al the available data set for the reaction it N — > tt tx N 
was analyzed within an isobar model, which allows for various intermediate 
7i 7i N states, such as A 71 and iV p. All partial waves were fitted with res- 
onance plus background terms, and from that fit the resonance parameters 
were extracted. Similar work has been done by various other groups jE] 



We follow the results of Manley et al, since his analysis is the most extensive 
and recent one. In this way, rather than by using the parameters of [pCfl , 
experimental information on pN scattering enters into the model in a well 
defined way. 

In Table [I] we list the mass uir, the total decay width T tot and the decay 
width into the pN channel, as taken from |19| . For comparison, we also give 
the value for the p N decay width from [^D| . In the last column we denote if 
the respective resonance decays in a relative s-wave or p-wave into the p N 
channel. 

A considerable uncertainty in the model is connected with the coupling 
strength of some resonances to the pN channel. This is most evident for 
transversely polarized p mesons. Here, at intermediate and high momenta, 
the Pi3(1720) resonance plays a dominant role, due to its large partial decay 
width into the p N channel as given in the analysis of Manley et al [ HJ . As 
listed in Table [l|, the estimate of the PDG group [p(J is smaller by a factor 
of 3. Moreover, in Manley's analysis a large error of 168 MeV - which is of 
about the same size as the decay width itself - is stated. Similarly, for the 
-D 33 (1700) the partial decay width as given by PDG and by Manley differ 
by a factor of 3, mainly affecting the results of A L p . Again, Manley gives a 
very large error of 31 MeV for this quantity. On the other hand, very good 
agreement is reached on the properties of the _Di3(1520). Here, Manley and 
the PDG group present nearly identical values for the partial decay width. 
Furthermore, the error for this number is much smaller, only 6 MeV are listed 
in@. 

In Manley's fit a few high lying resonances, which have not been included 
in our previous calculation, show a sizeable branching ratio into the pN 
channel. They are listed below the double-line in Table [TJ. However, as will 
be discussed in Sect. [5], they induce only small changes on the results. 

The total decay width V of the resonances is taken as a sum over all 
partial decay channels. The branching ratios into these channels as well as 



the mass dependence of the decays are the same as in [ 19| , except for the 



pN decay. This decay channel is calculated within the same framework as 



16 



the forward scattering amplitude, see Eq. |. 

Due to the small phase space there is no measured decay width of the 
P33(1232) into pN . Therefore the coupling constant of this resonance to the 
p N system can not be obtained by a fit to the decay width and is subject 
to substantial uncertainty. In our previous calculation, as well as in other 
works [[T6| , |17jl , it was taken from the Bonn potential model. It is however 
questionable, whether this value has much meaning outside the Bonn model 
(cf. discussion in [fill ). We therefore calculate the coupling constant directly 
from the well known isovector part of the A — > N'y decay width, assuming 



the validity of strict VMD |I2" |. Note that the numerical values obtained in 
both ways are nearly identical. The mass dependence of the decay width is 
obtained in the same way as for the other resonances. 

It might be argued, that the consistency of the resonance parameters 
employed for the calculation of the p N forward scattering amplitude can be 
tested by comparing with data on the reaction 7r iV — > p N |[43|| . On inspection 
of these data, there seems to be a lack of strength below the threshold at 
roughly 1.7 GeV, in conflict with a large coupling of the _Di3(1520). However, 
this conclusion is not correct. The pN "data points" were extracted from 
the measured n it N data by simple fits of the mass differential cross section. 
Naturally, with such a procedure no sensitivity is achieved in the subthreshold 
region, where the p meson can not directly be identified from the it it mass 
spectrum. Such a fit is certainly much less conclusive than the partial wave 
analysis of Manley et al. Moreover, the data set used by Manley et al also 



comprises the tt N — > tx tx N data from E3 



Results 1 - Relativistic Versus 
Non-Relativistic 



In Sect. 2.2 we discussed the shortcomings of the non-relativistic calculation 

T / L 

of the self energy S med of the p meson. We then developed the framework 
for a fully relativistic calculation. Now we compare the non-relativistic and 
the relativistic approach - using the relativistic Lagrangians from Eq. [17| - 
and pay special attention to the question, whether a reference frame exists, 
in which the non-relativistic reduction should be performed in order to arrive 
at a good approximation of the relativistic results. 
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In a resonance model, neglecting level repulsion for the moment, the 
spectral function A p is most sensitive to the behaviour of the scattering 
amplitude in the vicinity of s = m 2 R , i.e. along the lines depicted in Fig. 
|3|. In order to get a feeling for the differences between a relativistic and the 
different non-relativistic descriptions, it is therefore rewarding to compare 
the results from both approaches for the quantity Q T / L , which relates to the 
forward scattering amplitude as: 

T^ = - S/ F(A; 2 ) 2 - . (31) 

k — mi, + j v k 1 



Here k 2 stands for the invariant mass of the pN system; the formfactor F 
and the isospin factor Ej were already introduced in Sect. [|. Note that the 
coupling constants /rn p are included in the definition of Vt T l L and have been 
adjusted to the partial decay width of the resonance. The quantities Vt T / L 
are explictly given for both a non-relativistic and a relativistic calculation in 
the Appendix. 

In Fig. |] we show Q T for two of the most prominent resonances in the 
non-relativistic analysis, the -Di 3 (1520) and the Pi3(1720). It is plotted as 
a function of the invariant mass of the p meson, the momentum being fixed 
by the condition s = m 2 R . We compare the relativistic results with three 
non-relativistic versions of which are explained in detail in Sect. |T2 



i) relativistic 

ii) non-relativistic, selfenergy lab-frame, coupling constant cm-frame 

iii) non-relativistic, both quantities cm-frame 

iv) non-relativistic, both quantities lab-frame 

Up to now, in the literature only calculations performed within kinematics 



ii) can be found 0. 10 



We first remark that the results are not degenerate at threshold, since in 
each case the coupling constant has been adjusted to the measured partial 
decay width. As expected from the discussion in Sect. |2~2] , the calculations 
reveal large differences between the different non-relativistic approaches ii)- 
iv) in the transverse channel. This effect is - as already discussed - more 
pronounced for p-wave resonances, see lower graph of Fig. |1 Here, especially 
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at low invariant masses, the commonly used non-relativistic version ii) is not 
reliable at all, overestimating the results by a factor of 2 — 3, whereas both 
version hi) and iv) yield a surprisingly good agreement with the fully rela- 
tivistic calculation i). In the s-wave channel - upper graph of Fig. ^ - it is 
again the non-relativistic version ii) displaying a considerable disagreement 
with the relativistic calculation. Both versions hi) and iv) yield better ap- 
proximations. The exact form of the momentum dependence, however, is not 
well reproduced in version iv). We note that also in jPJ a good agreement 
between a relativistic calculation for the dilepton decay of the _Di3(1520) and 
a non-relativistic one performed in the cm-frame was reported. 

We do not show results for Q L . Here the effect from a relativistic cal- 
culation is rather small, since the momentum dependence of Q L is mainly 
determined from physical constraints at q = and at the photon point: at 
small momenta a relativistic calculation and its non-relativistic reduction 
must display a similar behaviour, whilst at the photon point Q L has to van- 
ish in either version. Thus, both calculations do not have much room to 
develop differently. There is a finite contribution to Q L from p-wave reso- 
nances, which, however, is small as compared to the transverse coupling. In 
the spin-i sector, one obtains qualitatively similar results for Vt T / L . 

Summarizing these results, we find that the previously used non-relativistic 
version ii) is not a good approximation of a fully relativistic calculation. In 
contrast, using cm- kinematics for both the determination of the coupling 
constant and the self energy leads to a very reasonable agreement with the 
relativistic approach. 

For this reason we argue also that the contribution from spin-| particles — 
for which a relativistic theory is very complicated — can be reliably evaluated 
in the non-relativistic approach hi). Since in a non-relativistic calculation the 
vertex factors of both the spin-| and spin-| resonances are proportional to 
q 2 , we expect from Fig. ^that at high momenta the contribution from spin-| 
resonances is reduced as compared to previous works [^, [TT], M, where q 2 
was evaluated in the lab-frame. 

Turning now to the discussion of the results for the spectral function, 
we start with a comparison between a relativistic calculation and the non- 
relativistic approach with kinematics hi). These and all other calculations 
of the spectral function have been performed at normal nuclear matter den- 
sity po = 0.17//m 3 . The good agreement for the results of f2 T//i translates 
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into very similar results for the spectralfunction Ay L , see Fig. [5| There the 
transverse spectralfunction is plotted at momenta q = GeV and q = 0.8 
GeV. The resonance parameters are taken from Manley et al, only the reso- 
nances above the double line in Table |] are included. The small deviations at 
large momenta are partly due to the different choices for the formfactor, see 
Eqs. |21] and |22|. The effect of the formfactor will be discussed in more detail 
later on in this Section. We do not present a comparison for the longitudinal 
channel as there the differences are even smaller. 

Since the agreement between these two approaches is nearly perfect we 
will use the remainder of this Section for the comparison of our previous 
non-relativistic results (within kinematics ii)) and the relativistic ones. Dif- 
ferences between both versions are not only due to kinematics. Also the 
changes in the resonance parameters and the influence of the new formfac- 
tor need to be considered. This is apparent already at low momenta, see 
Fig. |6]. For the moment we concentrate on a comparison of our previous 
non-relativistic calculation (dotted line) with the relativistic one, where the 
same resonances as in the non-relativistic case (listed above the double-line 
in Table [I]) are included (dashed line). At q = the non-relativistic calcu- 
lation with PDG parameters shows a stronger depletion of the p peak than 
the relativistic approach with Manley's parameters. 

The formfactor F(s) in Eq. ^1] suppresses off-shell contributions (s ^ m 2 R ) 
of resonances and thus ensures that at large invariant masses of the p meson 
the forward scattering amplitude becomes negligible. It therefore provides a 
conservative estimate of the resonance contribution. Such a behaviour is not 
guaranteed by using the old formfacor Eq. ^2] instead, which is equal to one 
at q = 0. With this formfactor the forward scattering amplitude increases 
with the invariant mass of the p, due to the energy dependence of the vertex 
function VL T I L . This affects in general both the tail of the resonance excitation 
and, due to level repulsion, the strength in the resonance peak. Thus, due to 
the formfactor, the impact of the excitation of the Di 3 (1520) in the region of 
the p peak is noticeably reduced. Since the formfactor has such an impact on 
the results, we investigated the sensitivity on the chosen value for the cutoff 
parameter A by varying A in the range of 1 — 1.5 GeV. The effects on the 
results for the spectral function were found to be small, they are at the order 
of 10%. This effect is enhanced by the small partial p N decay width of the 
-D33(1700) found in the Manley analysis as compared to PDG (ct. Table [l]). 
Its small coupling makes this resonance play a merely negligible role in the 
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relativistic calculation. 

At larger momenta both approaches lead to very different results for A T , 
see Fig. [7]. As a general tendency, in the new calculation much more strength 
is centered around the original p peak. This behaviour can be retraced to 
various sources. As follows from the properties of f2 T//L , for a given branch- 
ing ratio, a relativistic calculation produces a much smaller T^ ed than a 
non-relativistic one. The resonance contribution is further suppressed by 
the formfactor F(s), thus producing a more conservative estimate of the 
resonance contribution. In addition, at higher momenta, S med is reduced 
by the smaller branching ratio of the _D33(1700). These depletion mecha- 
nisms are only counteracted by the strongly increased branching ratio of the 
Pi3(1720). Overall, at high momenta A T p is dominated by the £)i3(1520) 
and the Pi 3 (1720) as well as, to some extent, by the spin-| resonances 
Pis (1680) and i<35(1905). Still the transverse spectral function receives sig- 
nificant broadening with increasing 3-momentum in a relativistic calculation, 
but this effect is not as pronounced as in the non-relativistic case. 

Whereas A T receives strong qualitative modifications in a relativistic cal- 
culation, changes in A L are only of a quantitative nature. This is shown 
in Fig. |8]. As already mentioned, the momentum dependence is dominated 
- both in a relativistic and a non-relativistic approach — by physical con- 
straints at q = and q 2 = 0. Consequently, the Di 3 (1520) peak at q = 0.4 
GeV contains nearly the same strength in both approaches. The depletion 
of resonance strength around the p peak is, as in the case of q = 0, due to 
the formfactor and the small coupling of the P 33 (1700). In the longitudinal 
channel the only sizeable contribution comes from the Z?x3(1520). 

The impact of the new, high lying resonances on our results is small and 
noticeable only at high momenta. This is easily understood for the following 
reasons: Despite their large branching ratio, the coupling strength is rela- 
tively small, since their mass is far above the p N threshold and consequently 
the phase space is very large. In addition, they have a very large total decay 
width, leading to a further suppression of the contribution. At low momenta, 
these resonances can only be excited from p mesons far above their pole mass. 
Correspondingly, their contribution to the imaginary part of £ me( 2 is drowned 
by the imaginary part of the vacuum self energy, the 2 7r decay width. It is 
only at momenta around 0.8 GeV, that the resonance excitation occurs at 
invariant masses around the pole mass of the p meson. At this point reso- 
nance scattering becomes important. This is shown in Figs. |6|-|8], where the 
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calculation including all resonances is represented by the solid line. 

Finally, we discuss the effect of medium modifications of the scattering 
amplitude. The Fermi averaging shows sizeable impact only at large mo- 
menta, leading to a broadening of the resonance peaks. As for the case of 
Fermi averaging, short range correlations influence the results mostly at high 
momenta. We show a comparison of the influence of both modifications on 
A T p at q = 0.8 GeV in Fig. p. For A L p the differences are completely negligible. 

6 Results 2 - Influence of Couplings 

So far, we have discussed the differences between a relativistic and a non- 
relativistic calculation of £ me <2 and their consequences for A p . However, the 
Lagrangians from Eq. |T7| is not unique; the resonance can couple in many 
more ways to the pN system. In order to draw safe conclusions on the 
momentum dependence of the pN scattering amplitude, it is important to 
study the influence of the chosen coupling on the results. 

In this analysis we confine ourselves to the case of spin-| resonances of 
both positive and negative parity. The vertex functions are derived from the 
Lagrangians in Eqs. pj] and [18|. In the spin-| sector there is no resonance 
which has a sizeable coupling to pN. Clearly, a strong dependence of the 
final results on the coupling of spin-| resonances is not expected. 

In the following discussion we will use some shorthand notations to dis- 
tinguish between the different vertex functions. The standard coupling from 
Eq. [I7|will be referred to as coupling i), to the first vertex function from Eq. 
|T8| we refer as coupling ii) and to the second one as coupling iii). To get a 
feeling for the changes, we present first VL T / L in Fig. ^ and |ll]. The results 
are shown for the _Di3(1520) and the Pi3(1720) resonances. Explicit results 
for these quantities are given in the Appendix. Note that due to dimensional 
reasons in the definition of VL T / L in Eq. [JT] an additional factor — has to 
be introduced, see all Eq. [18]. The results are nondegenerate at threshold 
since the coupling constants fnN P are already adjusted to the pN branching 
ratios. The most distinct behaviour is exhibited by coupling iii). At the pho- 
ton point, it enforces a vanishing contribution not only in the longitudinal, 
but also in the transverse channel. This behaviour can be read off from the 
Lagrangian £ int : 
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d v = -q 2 A v . (32) 

Clearly, the vertex function has to vanish at q 2 = 0, regardless of the 
polarization of the p meson, as can easily be read off the explicit expression 
for VL T / L in the Appendix, see Eq. [|2|. Consequently, S^ed nas t° fulfill the 
same constraints as S^ed- Both become more or less identical. Thus, the 
coupling iii) mostly affects the transverse channel. Noteworthy is also the 
much larger contribution from p-wave resonances to Q L . 

For the same p N branching ratios, large coupling constants are obtained 
by using coupling iii), as compared to the other two couplings schemes. In 
the determination of the coupling strength the quantity 2 Q T + Q L enters, 
averaged over the p mass distribution. Since this expression has to vanish 
at q 2 = in coupling iii), its size is much reduced at low invariant masses. 
This behaviour can only be compensated by a large coupling constant. For 
p-wave resonances this effect is most evident, since in this case the matrix 
element becomes sizeable mostly near the photon point for the couplings i) 
and ii). As an example, we obtain a value of fnN P = 36 for the Pi3(1720) in 
coupling iii) as compared to fnN P = 10 in coupling i). Clearly, the increased 
values for the coupling constants are responsible for the enhanced resonance 
contribution in the longitudinal channel. 

For coupling iii) the coupling strength for the P33(1232) can not be ob- 
tained from VMD since the matrix element vanishes at the photon point. We 
take an ad hoc value of fnN P = 30 for the coupling constant, which is twice 
as large as that obtained from coupling i). Varying fnN P between values of 
15 — 30 has only negligible influence on the results. 

Coupling ii) produces much less dramatic changes as compared to iii). 
The general properties of Vt T ^ L do not change much in comparison to a 
calculation employing coupling i). An exception is the coupling of p-wave 
resonances to longitudinal p mesons, where coupling ii) generates a much 
larger contribution than coupling i). Quantitatively, also the coupling of s- 
wave resonances to transverse p mesons change up to some extent. Nearly 
identical results are achieved in the remaining channels. 

Turning now to the results for the spectral function, at q = no signifi- 
cant modification is observed, see Fig. |T2j. However, at higher momenta, the 
distinct character of coupling iii) becomes apparent, reducing the resonance 
contribution at low invariant masses for transverse p mesons. At the same 
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time, a stronger modification of A L is exhibited. Thus, coupling iii) induces 
a qualitatively different behaviour on A p , removing much of the distinction 
between transverse and longitudinal p mesons. The results at high momenta 
are shown in Fig. |H| In this kinematical region spin-| resonances, for which 
we have only a non-relativistic coupling at hand, have noticeable influence 
on the results. Therefore, by changing the interaction in the spin-| sector 
and leaving the spin-| sector unchanged, the effect of the coupling scheme 
on the spectral function is underestimated. In order to obtain an upper limit 
for this effect, we treat the spin-| resonances the same way as spin-| reso- 
nances, correcting for the different spins of the states by a factor of 3/2. This 
procedure is reasonable, since — except for the multiplicity — the spin of 
the resonance affects mostly the angular dependence of the matrix element. 
We have checked that for the transition from spin-| to spin-| resonances our 
prescription works fine within 10%. 

The most distinct feature of coupling ii) is an increased resonance contri- 
bution to A L at large momenta. The transverse channel, however, does not 
undergo large modifications as compared to coupling i). 

One should keep in mind, that in general the actual coupling of a res- 
onance to the pN channel might contain contributions from all three cou- 
plings. Naturally this will smear out the differences. One expects from the 
above discussion, that the general tendency will go towards a less modified 
transverse p meson than predicted from coupling i), whereas in the longi- 
tudinal channel a stronger modification is likely. It is beyond the scope 
of this work, however, to find realistic weighting factors for each coupling. 
This can only be achieved by a complete relativistic partial wave analysis of 
7r N — ► 7T ii N data. 



7 Conclusions 

We have performed a calculation of the p spectral function A p in nuclear 
matter within a relativistic resonance model. The aim was, to put predic- 
tions concerning the properties of p mesons in nuclear matter, which were 
derived in a non-relativistic reduction ||10|| , on a more solid basis. Especially 
the momentum dependence of A p was examined. This discussion was mo- 
tivated by the observation that the results of a non-relativistic calculation 
depend strongly on the reference frame in which the non-relativistic vertex 
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functions are evaluated. We also investigated the influence of different res- 
onance parameters. We argued, that the resonance parameters as obtained 
in the work of Manley et al [IB], |TU[ provide much more reliable information 
on pN scattering than the estimated resonance parameters listed in PDG 
|p0|| . Finally, we studied if and to which extent the obtained results depend 
on the structure of the chosen coupling. Therefore, various coupling schemes 
for the RNp vertex were employed in the spin-| sector. 

A very interesting finding of this work is that a non-relativistic approach 
is in very good agreement with a relativistic one and should be good enough 
for all practical purposes, provided that the vertex functions for both the de- 
termination of the coupling constant and the scattering amplitude are eval- 
uated in the cm-frame, i.e. the rest frame of the intermediate resonance. 
In contrast, it is not advisable to determine the coupling constant with cm- 
kinematics but the scattering amplitude in the rest frame of the nucleon, as 
was done in all previous works || [10], [T^J. We could demonstrate that 



this procedure leads to dramatic deviations from a relativistic calculation. 

As an important finding in our previous non-relativistic calculation, we 
reported on the dominant role of the Di3(1520) in pN scattering. At low 
momenta these results are supported by the relativistic calculation. The 
parameters of this resonance are well known. Furthermore, different vertex 
functions show only marginal effects on the results. For higher momenta the 
role of the -D 13 (1520) is reduced in a relativistic calculation. 

Another striking outcome of a non-relativistic calculation was the very 
different behaviour of A T p and A L p at high momenta. In the light of a rela- 
tivistic version of A p this distinction needs to be revised. As an overall effect, 
taking into account all modifications on the calculations, a transversely po- 
larized p meson keeps a larger amount of its total strength around the pole 
mass than indicated in the earlier calculation. The previous finding of a com- 
pletely structureless strength distribution in the transverse channel at high 
momenta can therefore not be confirmed. In particular, the coupling with a 
derivative acting on the p field tensor produces a transverse spectral function 
A T , which undergoes less modifications at high momenta and is very similar 
to A^. 

The reliability of the calculation in the transverse channel is limited by 
poorly known resonance parameters, especially their partial decay width into 
the pN channel. In particular, a large uncertainty is connected with the 
branching ratio of the Pi3(1720). In this partial wave, experimental data do 
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not allow a conclusive determination of the resonance parameters. 

In the longitudinal channel only quantitative changes occur. Here, the 
results are governed by physical constraints at threshold and at the photon 
point. The change from a non-relativistic to a relativistic calculation has no 
substantial impact on the results. As in the transverse channel, at large mo- 
menta the coupling with a derivative acting on the p field tensor, produces 
the largest quantitative change. In this case it leads to a pronounced broad- 
ening of the p peak. Up to which extent this coupling is present in the actual 
RNp vertex remains an open question. As a general tendency, it will lead 
to a smaller modification of the transverse channel, whereas the influence of 
resonance scattering for longitudinal p mesons is enhanced. 

Altogether, we can confirm the dominant influence of the Di3(1520) on the 
propagation of p mesons in nuclear matter at low momenta. The sensitivity 
to the chosen coupling scheme is small and the physical parameters of this 
resonance are well under control. The distinction between A T and A L p as 
found in our non-relativistic study also exists in the relativistic calculation, 
though to a reduced extent. 
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A Appendix 



A.l \M 



RNp\ 



In this Appendix we show the problem arising from the completeness relation 
Eq. The squared matrix element l-M^^I 2 for the decay of a resonance 
with the quantum numbers J 71 " = | + into N p is readily calculated, using 
Eqs. [Li], [I7| and (erroneously) |26]: 



\M 



RNp\ 



P T 



{i>N + m N )Vi{^ + m R )Vl 
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4J r (~^~) mN (u 2 (mN + - m R ) + 



\M L RNp \ z = -P^Tr 



q (m N + m R -aj] 



(33) 




4 J r m N q (m N + uj - m R ) 



The calculation is done in the rest frame of the nucleon, and Jr denotes an 
isospin factor, see Eq. |2(| One easily verifies that \M 7 R ^ P \ 2 can become 
negative. At q = 0, for example, the squared matrix element changes sign 



at uj = m R — m N . Clearly, using Eq. [27] instead of Eq. g6| does not lead to 
such a pathological behaviour. 

For resonances with negative parity we do not find such effects. This 
can be understood by observing that the problem arises at the on-shell point 
y^s = m R . In the case of s-wave resonances a sign change at this point is not 
possible however, since the amplitude is no n- vanishing at the on-shell point. 
This, of course, does not imply that treating the s-wave resonances with the 
unmodified propagator is correct. 

A. 2 The Propagator 

According to |45| the propagator Di(k) has to satisfy a number of con- 
straints. Expressing the propagator in a Lehmann representation: 



Di(k) 



dM 2 



^p sp (M 2 )+M Psc (M 2 



k 2 - M 2 + ie 



(34) 



these constraints translate into properties of the spinor and scalar spectral 
weight functions of the propagator, p sp and p sc respectively: 



27 



i) p sp (M 2 ) and p sc (M 2 ) are both real 

ii) p sp (M 2 )>0 (35) 

iii) p sp (M 2 ) - Psc (M 2 ) > 

It is now possible to check whether the often used expression from Eq. 
|29| is compatible with these conditions. The comparison of Eqs. ^9| and |34| 
allows the identification of the functions p sp and p sc : 



i #r 



sp 



I rial 

1 m R T 



n(k 2 -m 2 R ) 2 + k 2 T 2 

p*: al (k 2 ) = ;;/;: w» 



n(k 2 -m 2 R ) 2 + k 2 T 2 
mR 



k 2 

Clearly p*" a ' and pf™ 1 , as obtained from Eq. ^ violate the last condition 
in Eq. for k 2 < m 2 



R- 

In order to construct a more reasonable propagator which incorporates 
the constraints from Eq. [35|, we will choose the spectral weight functions 
such that they fulfill these constraints and then calculate the full propagator 



directly from Eq. Keeping things as simple as possible we make the 



following ansatz for p sp and p sc : 

Ps P = Psc = P% al ■ (37) 

This is readily done for the imaginary part of the propagator, leading to 
the expression: 

ImDi(A0= Im ( ^ + \ _ (3g) 

5 \k 2 -m 2 R + iVk^Tj 

For the real part of Di (k), the integral can be solved analytically only, if a 
constant decay width F is assumed, leading to a lengthy expression. Taking 
into account the mass dependence of the decay width makes a numerical 
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computation mandatory. Concerning the calculation of the p meson self 
energy we checked by explicit comparison that the ansatz 



Re Di(k) — Re ( ) (39) 

produces a satisfactory agreement with the exact numerical solution. Small 
deviations can be seen, but these have virtually no effect on the spectral 
function of the p meson. Obviously, Eqs. [38| and ^ lead back to our guess 
Eq. |§. 

One can write down a Lehmann representation for the propagator of spin- 
| particles D^{p) in a straightforward extension of Eq. [M[ This expression 

2 

contains in principle eight spectral weight functions, four multiplying the 
spinor part of D^{p) and four its scalar part. By analogy to the case of 

2 

spin- 1 fields, we require them to be identical and put them equal to p 1 ™ 0,1 



from Eq. KB. Thus, we obtain the result from Eq. 30 



A.3 Results for Q T / L 

In this Appendix we give explicit expressions for the quantities VL T I L which 
have been introduced in Sect. [5]. They are evaluated in the rest frame of 
nuclear matter, thus u and q refer to energy and momentum of the p meson 
in the lab-frame and q 2 denotes its squared invariant mass. The invariant 
energy of the intermediate resonance y^s is given by ^fs = q 2 + m 2 N + 2m nuj. 
Note that in the spin-| sector singular terms ~ - show up, in accordance with 
the discussion in Appendix 2. Note that all results need to be multiplied by a 

factor -^r-. We first list the results obtained from the relativistic Laerangians 
in Eq. O: 
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denoted by ii) in Sect. |^ one derives: 
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The coupling denoted by iii) in Sect. | leads to: 



[} T = | g 4 (m^v + oj — y/s) 



Q 1 



n 1 



m N q 4 (mjv + u + y/s) ^1 
m N g 4 (mjy + ^ + 



m N \ q 2 

9 2 



(42) 



For comparison we also give the non-relativistic results | ID| . Note the 
appearance of an additional factor y/s + which needs to be introduced 
since in Eq. |ST| the scattering amplitude is formulated using a relativistic 
propagator. 
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(43) 
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Here q refers to either the cm-momentum or the lab-momentum of the p 
meson, depending on which kinematic has been chosen. 
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Table 1: List of the resonances included in this calculation. The parameters 
in the first three columns are taken from the Manley analysis |18[ [19]. In the 
fourth column the resonance decay width into pN as presented in the PDG 
[ pCfl is listed. The fifth column shows the coupling strength of the resonances 
to pN as obtained from expression @. In the last column we denote if 
the respective resonance decays in a relative s-wave or p-wave into the p N 
channel. The resonances below the double line have not been included in our 



previous analysis 10 




Figure 1: Resonance contribution to the pN scattering amplitude. 
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Figure 2: Decay of a resonance into two pions via an intermediate p meson. 
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Figure 3: Free dispersion relations for the p and the most prominent reso- 
nances. For the resonance branches m is the invariant mass a meson of given 
three-momentum q must have in order to excite the respective resonance on 



its mass shell |IOfl. 
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Figure 4: Q T for the D 13 (1520) resonance (upper graph) and the Pi 3 (1720) 
(lower graph) at s function of the invariant mass m of the p 

meson. Shown are three different calculations, which are explained in Sect. 

B 
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Figure 5: The transverse spectral function A T p at GeV and 0.8 GeV for 
p = po as a function of the invariant mass m of the p meson. Shown is a 
comparison between a relativistic calculation and and a non-relativistic one 
using cm-kinematics iii). 
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Figure 6: A T p = A L p at q = for p = po as a function of the invariant mass 
m of the p meson. Compared are a relativistic (solid line, dashed line) and 
a non- relativistic calculation (dotted line). For the relativistic calculation 
the resonance parameters of Manley et al |19[ are employed, for the non- 
relativistic one the parameters of PDG [20|. 
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Figure 7: A T p at q = 0.4 GeV (upper part) and at q = 0.8 GeV (lower part) 
at p = pO as a function of the invariant mass m of the p meson. 
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Figure 8: Same as Fig. [7|, but for A L p . 
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Figure 9: A T p at q = 0.8 GeV for p = po as a function of the invariant 
mass m of the p meson. The vacuum result (solid line) is compared to an 
in-medium calculation where Fermi averaging and short range correlations 
are taken into account (dashed line). 
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Figure 10: Q T (upper graph) and Q L (lower graph) for the .Di3(1520) at 
function of the invariant mass m of the p meson. Shown is 
the influence of different coupling schemes on the results. The couplings are 
explained in Sect. 
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Figure 11: Same as Fig. [TO but for the Pi3(1720) 
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Figure 12: A T p = A L p at q = GeV for p = p as a function of the invariant 
mass m of the p meson. The couplings employed for the calculations are 
explained in Sect. ||. 
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Figure 13: A T p (upper graph) and A L p (lower graph) at q = 0.8 GeV for 
p = po as a function of the invariant mass m of the p meson. The couplings 
employed for the calculations are explained in Sect. |6|. 
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